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ABSTRACT 

Magnetic fields, which are undoubtedly present in extragalactic jets and responsi- 
ble for the observed synchrotron radiation, can affect the morphology and dynamics 
of the jets and their interaction with the ambient cluster medium. We examine the 
jet propagation, morphology and magnetic field structure for a wide range of den- 
sity contrasts, using a globally consistent setup for both the jet interaction and the 
Q ■ magnetic field. The MHD code NIRVANA is used to evolve the simulation, using the 

constrained-transport method. The density contrasts are varied between rj = 10~^ 
C/3 , and 10~^ with constant sonic Mach number 6. The jets are supermagnetosonic and 

simulated bipolarly due to the low jet densities and their strong backfiows. The helical 
magnetic field is largely confined to the jet, leaving the ambient medium nonmagnetic. 
We find magnetic fields with plasma /? ~ 10 already stabilize and widen the jet head. 
^ ' Furthermore they are efficiently amplified by a shearing mechanism in the jet head 

lO \ and are strong enough to damp Kelvin-Helmholtz instabilities of the contact disconti- 

•/^ I nuity. The cocoon magnetic fields are found to be stronger than expected from simple 

Op ' fiux conservation and capable to produce smoother lobes, as found observationally. 

, The bow shocks and jet lengths evolve self-similarly. The radio cocoon aspect ratios 

00 ' are generally higher for heavier jets and grow only slowly (roughly self-similar) while 

' overpressured, but much faster when they approach pressure balance with the ambi- 

0^ . ent medium. In this regime, self-similar models can no longer be applied. Bow shocks 

' are found to be of low excentricity for very light jets and have low Mach numbers. 

Cocoon turbulence and a dissolving bow shock create and excite waves and ripples in 
the ambient gas. Thermalization is found to be very efficient for low jet densities. 

^ ■ Key words: galaxies: jets - magnetic fields - MHD - methods: numerical - radio 

' continuum: galaxies - galaxies: clusters: general 



1 INTRODUCTION 

Extragalactic jets are amongst the most powerful phenom- 
ena in the universe and observable up to high redshift 
l|Milev fc De BreuckllioO^ ). Jet activity is generally accom- 
panied by synchrotron emission from relativistic electrons in 
the magnetized jet plasma, which is most prominently ob- 
servable at radio frequencies. Thus radio observations pro- 
vided us with much insight into jet morphology (beams, 
knots, hotspots, lobes/cocoons), classifi cation into low- 
powe r FR I and high-power FR II sources (|Fanaroff fc RilevI 
1 1974 ). magnetic field strengths estimated by minimum en- 
ergy arguments in the ra nge of a few to several hundreds 
of microgauss in hotspots l|Bridlelll982l : iMeisenheimer et al.l 
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as well a s magnetic field topo logy from polarization 
measurements (|Bridle fc Perlevlll984l ) (high power jets gen- 
erally show magnetic fields paral lel to the jet axis) and age 
estimates from sp ectral ageing l|Alexander fc Leahvl Il987l : 
ICariUi et al.llT991l ) of a few x 10^ years. Extragalactic jets 
show prominent cocoons, often only partly visible as lobes, 
with total length to total wid th aspect ratios mostly between 
3 and 10 jMullin et al.ll2008h . However, it must be cautioned 
that visible radio emission does not necessarily trace all the 
bulk plasma, but depends on magnetic field strength and 
electron acceleration. 

The current X-ray observatories Chandra and XMM- 
Newton give rise to a complementary view on jets, ob- 
serv ing bremsstrahlun g radiation from the thermal ambient 
gas (ISmith et al ll2002l) and inverse-Co mpton emission from 
the jet cocoon IjHardcastle fc Crostonllioosi : [ Croston et al.l 
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I2005I ). the latter suggesting near-equipartition magnetic 
fields. These emission processes are less dependent on micro- 
physics and thus easier to connect with theoretical results 
and furthermore contain information about the history of 
these sources. 

An extreme exam ple for this is MS0735. 6+7421 
l|McNamara et ahllioosl ), where radio emission shows a weak 
source, while the spatially coincident X-ray cavities reveal 
the true average power of the AGN (1.7 x lO^^ergs"^), 
which is a factor of ~10^ higher. The jet cocoon displaces 
the ambient thermal gas and drives a bow shock outward. 
Both pV work and the bow shock contain signatures of 
the jet properties and may be excellent diagnostic tools. 
Nearly three dozen cluster s were found to have cavities 
l|McNamara fc Nulsen|[2007l ) and generally show weak bow 
shocks (Mach 1-2) with aspect ratios (length/width) not 
much above unity. 

Radi o observations of h igh-power FR II jets show wide 
cocoons ijMuUin et al.l [2008I ). which are partly (and more 
completely at low frequencies) visible as radio lobes. Both 
jet head and cocoon-ambient g as interface appea r smoother 
in radio maps (e.g. Cygnus A (jLazio et al.ll2006l'). Pictor A 
llPerlev et al.lll997l ). but also Hercules A l|Gizani fc Leahvl 
[2003)) than they do in hydrodynamic simulations; and while 
the emission of synchrotron radiation obviously indicates the 
presence of magnetic fields, their importance for jet dynam- 
ics and the contact surface is still unclear. 

On the the oretical side, early n umerical simulations of 
supersonic jets l|Norman et al.|[l982l ) already exhibited basic 
structures seen in observations of extragalactic radio sources 
(working surface, cocoon, bow shock) and showed that pro- 
nounced cocoons are properties of jets with much lower den- 
sity than the ambient medium (light jets), although the 
slow propagation of the jet makes simulations of these com- 
putationally very expensive. With the availability of more 
computing pow er and new codes, sim ulations of the long- 



term evolution jReynolds et al. 20021) . in three dimensions 
llBalsara fc Norman! 1 19921 : IClarke et al]|l997j : lKraus3l2005l: 



Heinz et 3.1.1120061) ■ of very hght (lKrausel l2003': 'Saxt on et all 
'2002a, b:~ICarvalho fc O'Deal l2002al lbl; Izanni c t al.l l2003l: 
Sutherland & Bickncll 2007 ) and relativistic jets (|Alov et al.l 
1999»: ,Komissarovi ,199a : .Rosen et al.l 1 19991 : iHarded \20^ . 



Leismann et al.l [2005|) became possible . Furthermore, ef- 



fects of magneti c fields were examined (Clarke ct al. 1986; 
Lind et al.lll989l: iKossl et al.|[l990: Hardee fc Clarke 1995 : 



Tregillis et al. 



Mizuno et al. 



2001 



20041: lO'Neill eit al.ll2005 l: iLi et aljbood : 



20071 : iKeppens et al.ll2008l ~ although only for 
relatively dense jets (density contrast ^ 10~^). In this paper, 
we extend the studies of very light jets to include magnetic 
fields. 

To explore the interaction of jets with the ambient intra- 
cluster medium and and the impact of magnetic fields, we 
performed a series of magnetohydrodynamic (MHD) simu- 
lations of very light jets on the scale of up to 200 kpc (200 
jet radii) with a globally consistent magnetic field configura- 
tion. A constant ambient density was used to avoid effects of 
a declining cluster gas density on the structural properties 
which could contaminate the effect of magnetic fields, while 
the effe c ts of a density profile previously were described in 
iKrausd (|2005l ) for the axisymmetric and threedimensional 
case. 

After a description of our simulation setup and the nu- 



merical method, our results are described: first about the 
morphology and dynamics, the evolution of the bow shock 
and the cocoon; then entrainment of ambient gas and the 
energy budget; and finally the magnetic fields and their evo- 
lution as well as their impact on morphology and propaga- 
tion. Results are then discussed and put in context with 
observational findings. 

To compensate for different propagation speeds of jets 
with differing density contrasts, plots will use the axial 
bow shock diameter or jet length, where appropriate. "Full 
length" refers to the whole simulated length (considering 
both jets), while "full width" refers to twice the measured 
(radial) distance from the axis. 



2 SETUP AND NUMERICAL METHOD 
2.1 Setup 

The idea behind the present study was to explore the be- 
haviour of very light jets with non-dominant magnetic fields 
in a cluster environment using a plausible global setup for 
the plasma and the magnetic fields, but still keeping the 
setup simple enough to see the working physical processes 
clearly, which is much harder for a complex setup. We 
performed 2.5D simulations (axisymmetric with 3D vector 
fields) of both purely hydrodynamic and MHD jets on the 
scale of up to 200 kpc with a constant ambient gas density, 
where density contrasts -q = pj/pa were varied between 10~^ 
and 10""* to see its effects on the simulation. Jet speed, beam 
radius, sonic Mach number and ma gnitude of the helical 
magnetic field (|Gabuzda et ahlboOSl ) were kept fixed, thus 
yielding a kinetic jet power Lkin = nr^rfp^v^ and plasma 
P = ^■Kp/B'^ varying with density contrast. A summary of 
the parameters is given in Tab. [T] and [21 The simulations 
are labelled by a letter and a numeral, indicating the inclu- 
sion of magnetic fields (M) vs. pure hydrodynamics (H) as 
well as their density contrast. Both simplifications - axisym- 
metry and density di stribution - w ere relaxed in a previous 
hydrodynamic study (jKraus and their infiuence is ad- 

dressed later in the discussion. The initial gas distribution 
was randomly perturbed on the resolution scale to break 
symmetry between both jets. 

The bipolar (back-to-back) jets were injected by a cylin- 
drical nozzle (radius rj , length 2rj ) along the Z axis in cylin- 
drical coordinates {Z,R,(j)), hence allowing for interaction 
of the backflows in the midplane. The jet radius is resolved 
by 20 cells. Fully ionized hydrogen (7 — 5/3) was assumed 
for both the jet and the external medium. A compressible 
tracer field was advected with the fiow, using a value of 1 
for the ambient gas, and and —1 for the jets, allowing to 
trace back the origin of the plasma. Optically thin cooling is 
included in the code but was switched off, since the cooling 
times for our setup (S> 10* years) are significantly longer 
than the simulated time-scale, even for the shocked ambient 
gas. This choice also makes the simulations scalable, e.g. to 
other values of the jet radius, which was chosen arbitrarily 
as rj = 1 kpc. 

In the jet nozzle, all hydrodynamic variables (density, 
pressure and velocity) are kept constant at all times and 
a toroidal field oc sgn(Z) sin'*(7ri?/rj) i?/rj is prescribed 
there, being zero outside the nozzle. A dipolar field cen- 
tred on the origin is used as initial condition for the whole 
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Table 1. Common simulation setup parameters. Nozzle averages 
are restricted to _R ^ 0.9 rj 0.8 rj for M4 and M4L). 



jet speed 




0.6 c 


jet sound speed 


Cs 


0.1c 


jet radius 




1 kpc 


ambient gas density 


Pa 


IQ-^ nip cm~^ 


ambient gas temperature 


Ta 


5 X lO^K 


jet nozzle magnetic field 


(Bp) 


18.1 lO-G (M4L: 1.81 yiG) 






7.5 ^iG (M4L: 0.75 yiG) 



Table 2. Parameters for the difi^erent simulation runs. Nozzle 
averages are restricted to i? ^ 0-9 ''j 0.8 rj for M4 and M4L). 
/3 values are typical over simulation run. 



Run 


ri = Pj/pa 




tmax [Myr] 


Ml 


10-1 


810. 


6.7 


M2 


10-2 


81. 


10.9 


M3 


10-3 


8.1 


16.5 


M4 


lo-'' 


0.89 


47.5 


M4L 


10-4 


36. 


50.0 



computational domain (magnetic moment aligned with the 
jet axis) although it is mostly confined to the jet due to 
the strong decrease in magnitude with distance. For global 
simulations, the V • B constraint enforces closed field lines, 
which is satisfied by a dipolar field configuration, but not 
by the common setup of an infinite axial field, which lo- 
cally, but not globally fulfills the constraint. Thus, in our 
setup, the magnetized jet plasma propagates into the essen- 
tially nonmagnetic ambient matter. For M3 and especially 
M4, the magnetic fields become dynamically important and 
infiuence the appearance, so another run with lowered mag- 
netic fields was performed in addition (M4L) , which is more 
in line with the other jets. These lightest jets are addressed 
more specifically in Sect. 13.71 

The initial magnetic fields in Tab.[T]are nozzle-averaged 
initial values. As the poloidal field cannot be kept constant 
in the nozzle without violating V ■ B — 0, this field can 
evolve with time due to the interaction with the enclosing 
cocoon, quickly adjusting to 13 [iG (1.5 ^.G for M4L) but 
then stays constant. For these nozzle averages, only 90 per 
cent of the jet radius were considered for Ml, M2, M3 and 
80 per cent for M4 and M4L to exclude cells at the shearing 
boundary of the nozzle, where high magnetic fields and op- 
posite field directions can occur, while the core of the jet is 
unchanged. For plasma (3, we use the volume- averaged har- 
monic mean, as very weakly magnetized regions otherwise 
would misguidingly dominate the average. 

Since the jets are very underdense with respect to the 
ambient gas and have a high internal sound speed, a re- 
confinement shock develops already very near the jet noz- 
zle. This shock establishes pressure balance between the jet 
beam and the cocoon, resulting in a pressure-confined beam. 
In contrast to freely expanding jets, any imposed opening 
angle becomes unimportant once the beam is reconfined al- 
ready near the jet inlet. Hence, in contrast to heavier jets, 
underdense jets quickly find pressure balance with their en- 
vironment. 

The simulations were run until they reach the boundary 



of the uniform grid which has 4000 x 800 or 4000 x 1600 cells 
(for M4 and M4L) and the jet radius (rj) is resolved with 
20 cells. 

In the following, we will focus only on the MHD jets, 
as their hydro counterparts are only for comparison (set up 
exactly as the MHD jets with vanishing magnetic field). 



2.2 Method 

The simulations were p erformed using the NIRVANA code 
IZiegler fc Yorkelll997l ). which numerically solves the non- 
relativistic magnetohydrodynamic equations in three dimen- 
sions in cartesian, cylindrical or spherical coordinates. It is 
based on a finite-differences discretization in an explicit for- 
mulation using operator splitting and uses van Leer's in- 
terpolation scheme, which is second order accurate. The 
advection part is solved in a conservative form and the 
magnetic fields are evolved using the constrained transport 
method, which conserves V ■ B up to machine roundoff er- 
rors. T he code was vector ized and shared-memory paral- 
lehzed (jGaibler et al.ll2006l ) for the NEC SX-6 and SX-8. 



3 RESULTS 
3.1 Morphology 

The density and temperature maps of Fig.[T]show snapshots 
of all runs at a full jet length of 100 kpc, respectively. In the 
following, M3 mostly will be used for figures as it has the 
strongest non-dominant magnetic fields, therefore showing 
effects of the magnetic fields best and allowing for compari- 
son of features between different figures. 

The jet backflow blows up a pronounced cocoon, sur- 
rounded by a thick shell of shocked ambient matter. Dense 
ambient gas is mixed into the cocoon in finger-like struc- 
tures due to Kelvin-Helmholtz instabilities at the contact 
surface. Near the jet heads, this instability is suppressed by 
the magnetic field, which leads to a smoother appearance 
there. In purely hydrodynamic simulations, this stabiliza- 
tion is absent. 

The cocoon is highly turbulent and vortices hitting the 
jet beam can easily destabilize, defied or disrupt it if jet den- 
sities are l ow. The Mach number varies consi derably along 
the beam (|Krause|[200^ : ISaxton et al. l l2002bD and there is 
no stable "Mach disc" as seen for heavier jets - the terminal 
shock moves back and forth and often is not clearly defined. 

Because very light jets only propagate slowly, basically 
hitting the ambient g "solid wall", the backfiow is 

strong and the turbulence makes the interaction between 
both jets in the midplane important. Such jets have to be 
simulated bipolarly to describe the lateral expansion and 
hence the global appearance correctly. If only one jet were 
simulated, for very light jets the result would strongly de- 
pend on the bound ary condition in the equatorial plane 
l|Saxton et al.ll2002bl ). 

The surrounding ambient gas is pushed outwards by the 
cocoon pressure, driving a bow shock outwards. The bow 
shock for very light jets is different in its shape and strength 
from that of heavier jets (see Sect. I3.3|l . It is addi tionally 
changed by a density profile in the external medium (jKrausd 
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Figure 1. Density and temperature for all runs, eacli at a jet length of falOO kpc. The upper halves of each panel show the logarithmic 
density in units of 10~'^* g/cm^, the lower halves show the logarithm of temperature in units of 10^" K. Panels additionally are labelled 
by their respective time. 



[2OOI), which increases the aspect ratio with time because 77 
increases at the jet head and thus shows cyhndrical cocoons. 



3.2 Defining the cocoon 

In the following, we not only measure properties of the bow 
shock, which is easy to pin down, but also of the cocoon. 
While generally we define the cocoon as the region, which 



is filled by jet-originated matter (not including the beam 
itself) this definition has to be made in more detail for 
the simulation analysis. The strong backflow and the fragile 
beams of very light jets make the distinction between cocoon 
and beam difficult, while mixing at the contact discontinu- 
ity complicates the assignment of cell to cocoon or ambient 
matter. While we do not attempt to distinguish between 
beam and cocoon if not stated explicitly (it only seems nec- 
essary for energetic investigations), the distinction between 
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Figure 2. Cocoon mass depending on tracer field threshold (mass 
of all cells with tracer below the threshold), for simulation M3 at 
t = 15 Myr. 



2.5x10"^ 




toroidal field limit [G] 

Figure 3. Cocoon mass depending on toroidal field threshold 
(mass of all cells with toroidal field magnitude above the thresh- 
old). Simulation M3 at t = 15 Myr. 

cocoon and ambient matter is necessary especially for the 
entrainment measurements later and thus is described in 
more detail in this subsection, along with the measurement 
of the cocoon properties. 

3.2.1 Cell assignment 

Two properties can be used for the distinction between co- 
coon and ambient matter: the (compressible) tracer field and 
the toroidal magnetic field. Tracer field values of 1 and above 
indicate undisturbed and shocked ambient matter. This is 
available in all simulations, but mixing with jet matter at 
the contact discontinuity (due to finite resolution) lowers the 
tracer and thus requires a threshold value. The cocoon mass 
is especially sensitive to this threshold value, as the den- 
sity of the ambient gas is much higher and thus causes large 
changes of the cocoon mass if the border is shifted. Figure[2] 
shows the cocoon mass for a range of tracer thresholds. The 
injected mass at this time is only 4x10^ Mq ; measured mass 
above this value is the entrained ambient gas mass. 

In contrast, the toroidal field strength can be used for 
separation, as the toroidal field is zero initially in the am- 
bient medium and is conserved independent from the other 
field components. Figure[3]shows the cocoon mass depending 




10 

time [Myr] 



Figure 4. Evolution of the cocoon width for different width def- 
initions for simulation M3. A tracer limit of 0.5 was used for cell 
assignment. 

on the toroidal magnetic field threshold. Using this method, 
even cells with only a small mass fraction of jet matter can 
be assigned to the cocoon. There is a clear break visible, but 
the cocoon mass continuously increases for lower threshold 
values until machine accuracy is reached. This has two ma- 
jor problems: First, it naturally is not available for the pure 
hydro simulations and thus cannot be used to compare HD 
with MHD simulations. Second, the high sensitivity to jet 
matter is not a real advantage, as the mass values do not 
nicely converge and we have to choose a threshold. 

In the following we will use a tracer threshold of 0.5 
which is available for HD and MHD models and gives cocoon 
masses that do not strongly depend on the tracer threshold. 
Furthermore, it selects the regions one would consider co- 
coon also by looking at the other physical variables. 

3.2.2 Shape measurement 

To characterize the width of the cocoon, we checked four 
different measures, which will generally give different results 
due to the ragged shape of the contact surface. Widths are 
measured from the symmetry axis {R = 0, jet channel) and 
thus are only "half widths" . Figure U shows the temporal 
evolution of the cocoon width definitions 

(i) maximum width: measured at the maximum R posi- 
tion of a cocoon cell, 

(ii) average width: Z-averaged over the full jet length, 

(iii) QB width: measured at one quarter the full jet length 
backwards of the jet heads, 

(iv) spheroid width: semi-minor axis of a spheroid with a 
volume equal to the cocoon volume and the semi-major axis 
equal to half the full jet length. 

The QB width is clearly much dependent on vortices 
near the contact surface and not very straight. Despite that, 
it grows similar to the Z-averaged width, which mostly has 
the lowest width value of all four definitions. The spheroid 
width lies between the maximum width and the Z-averaged 
width. All of theses measures can be approximated by power 
laws, although with somewhat different parameters. 

In contrast to the cocoon mass, the cocoon shape does 
not depend strongly on the tracer limit that is used for its 
determination (Fig. [5]). For limits around 0.5, the differences 
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Figure 5. Cocoon width for different width definitions as func- 
tion of the used tracer limit. For simulation M3 at 15 Myr. 
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Figure 6. Pressure maps in logarithmic scaling for simulation 
Ml and M3 at equal lengths. 



between width definitions is larger than the dependence on 
the tracer limit. 



3.3 Evolution of bow shock and cocoon 

3.3.1 Cocoon pressure evolution 

The low jet density has two main consequences for the evo- 
lution of the cocoon pressure: one is the lower jet power 
(for a fixed jet bulk velocity), which results in a lower co- 
coon pressure and a generally weaker bow shock. The other 
is the slow jet head propagation, which lowers the propa- 
gation time-scale compared to the dynamical time-scale of 
travelling pressure waves within the cocoon. Pressure waves 
from the jet head together with waves induced by turbulent 
motion and mixing in the cocoon, try to establish pressure 
balance within the cocoon and between cocoon and ambient 
gas, driving the lateral expansion of the cocoon. 

Figure |6] shows pressure maps of jets with r] — 10~^ and 
rj = 10~^ at the same lengths. The cocoon of the heavier jet 
is overpressured by a factor of 20 with respect to the ambient 
gas, while being a factor of only 1.5 for the hghter jet (and 
4.9 for this jet at the time of the Ml image). 

The strong evolution towards pressure balance is re- 
sponsible for the much less pronounced high-pressure re- 




-8.0 



E 

0) -9.0 



10.0 



M3 ® t 



-28 -27 -26 
log density [g cm"^] 

: 1 5 Myr 




66 



63^ 



-30 



-29 -28 -27 -26 
log density [g cm"^] 



-25 



Figure 7. Pressure vs. density histogram for the M3 jet after 2 
and 15 Myr. Volume-weighted counts, pressure and density are 
shown logarithmically. 



gions between Mach disc and the advancing bow shock. The 
bow shock has an elliptical shape with less directional de- 
pendence of its strength, more similar to an overpressured 
bubble, although it is still stronger in axial direction (see 
Sect. 13X211 . 

The quick pressure adjustment can also be seen in the 
pressure-density diagrams of Fig. [T] The ambient gas is de- 
scribed by the patch near (—26, —10), the jet nozzle by the 
cells around (—29, —10). Adiabatic compression and expan- 
sion leads to the oblique and longish features present at dif- 
ferent positions. Top right of the jet nozzle position are the 
cocoon grid points, which spread over a large range of den- 
sity to the right because of mixing with shocked ambient gas, 
which is the elongated feature top right of the ambient gas 
position. Comparing the two different simulation snapshots, 
we find that the pressure distribution is quickly adjusting to- 
wards the ext ernal pressure, in agreement to the findings in 
iKrausd l|2003l ). and the cocoon is not strongly overpressured 
anymore. 

Another view on this is the average cocoon pressure, 
shown in Fig. [S] which has a power-law-like behaviour. For 
the three models Ml, M2 and M3, it strikingly decreases 
with the reciprocal jet length (Tab. [3]). While Ml at the 
end of the simulation is still very overpressured, the cocoon 
pressure of M3 is already near the ambient gas pressure. 

M4 and M4L seem to deviate from this behaviour. At 
the beginning this is mainly a consequence of the longer- 
lasting relaxation from initial conditions, where strong 
shocks are thermalized efficiently, increasing the pressure in 
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Figure 9. Linearly scaled pressure map of M3 at 15 Myr. Values 
above 5 X 10"^" dyne cm~^ are clipped. 
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Figure 8. Evolution of the volume-averaged cocoon pressure as 
function of time or the jet length for the different models. The co- 
coon is defined by a tracer limit of 0.5. MHD models are shown in 
thick lines, corresponding hydro models in thin lines, the ambient 
pressure is indicated by the long-dashed line for comparison. 



the early "cocoon bubble" and because the early phase is 
shown with higher time resolution. After a jet length of 20 
kpc has been reached, they fit into the behaviour of the other 
simulations, but, as they soon reach the ambient pressure, 
settle to its value. 

It is clear that the cocoon pressure cannot drop much 
below the ambient pressure and thus approaches its value. 
At this point we expect the bow shock to softly turn into 
an ordinary sound wave. This is just about to happen in the 
last snapshots of M4 and M4L, where there is only a very 
weak density jump, corresponding to Mach 1.05. The exact 
value of the average cocoon pressure is insensitive to the 
exact definition of the cocoon (see Sect. [5^ . but can drop 
slightly below the ambient pressure due to pressure variation 
within the cocoon (which can still be as strong as a factor 
of 2). 

The past bow shock is not the only sound wave tes- 
tifying to the expanding cocoon. Already much before the 
shock decays, waves and ripples can be seen in the shocked 
ambient gas (Fig. [5J . 



3.3.2 Bow Shock 

The quick decrease in cocoon pressure naturally affects the 
strength of the bow shock as it is this pressure that drives the 
shock laterally. Figure [10] shows the temporal evolution of 
the bow shock strength, in terms of external Mach numbers. 




100 



axiol bow shock radius [kpc] 



10 S 



M3 ■ 
M4L ■ 




Qxiol bow shock radius [kpc] 

Figure 10. Evolution of the forward (top panel) and sideways 
(bottom panel) bow shock strength as a function of the mono- 
tonically increasing axial bow shock radius. Thick lines are MHD 
models, thin lines the corresponding hydro models for compari- 
son. 



for the forward direction (at R = 0) as well as the lateral 
direction (at Z — 0) for jets with different density contrasts. 

The bow shocks in forward direction are always stronger 
than the sideways shocks due to the direct impact of the jet 
on to the ambient gas. The lighter jets have a much weaker 
bow shock in all directions and the differences between the 
axial and lateral direction are much less pronounced. 

The axial diameter of the bow shock grows as a power 
law with exponents ~ 0.65 fFig. llll and Table[3]). For the lat- 
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Table 3. Power law exponents for the bow shock and cocoon evolution. The exponents p in the table are fits to power laws as function 
of time (for the lengths and widths) or as function of full jet length (for the average cocoon pressure), considering only the data points 
which match the criterion in the third line. A minus (-) denotes that no reasonable power law fit could be done. 
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Figure 11. Propagation of the bow shock in axial (top panel) 
and lateral direction (bottom panel) for the simulated models 
as function of time. Thick lines are MHD models, thin lines the 
corresponding hydrodynamic models. A power law oc t" '' is shown 
for comparison as long-dashed line. 



eral propagation we find similar exp onents. Th is behaviour 
agrees with self-similar jet models (iFalld Il991: Bcg clman 
1 19961 : iKaiser fc Alexandeij 1 19971 : iKomissarov fc FaUel Il998l ) 
and the spherical blastwave approximation ( Krausell2003l ). 
which predict an exponent of 0.6. At very early times and 
lasting longer for the lighter jets, we find lower exponents, as 



for a Sedov blast wave {I oc t ) from the initial conditions. 

The lighter jets have generally lower Mach numbers, as 
their kinetic power is lower, thus showing smaller bow shock 
velocities. The aforementioned analytical models yield an 
expansion speed 



TV Pa 



1/3 



-2/3 



kVjTI 



1/3 



-2/3 



(1) 



that directly translates into the bow shock Mach number 
and describes the scaling behaviour of the simulations rea- 
sonably well (Lj: jet power). We find values for k between 
1.5 and 2 in axial direction and between 0.5 and 2 (M4L) 
in lateral direction, the latter being increasingly higher for 
lighter jets. 

A clear deviation from this behaviour is M4 at t > 10 
Myr in axial direction. The bow shock propagates much 
faster due to the formation of a nose cone. The strong 
toroidal magnetic field coUimates the jet, suppresses the pro- 
nounced backflow of M4L, and the Lorentz force of the radial 
current gives the jet additional thrust for the propagation 
(see l3.7p . The other light jets (M3 and M4L) also may prop- 
agate somewhat faster due to their appreciable magnetic 
fields. 



3.3.3 Cocoon 

As the jet pushes the bow shock forward in axial direction, 
the cocoon length grows similar to the bow shock length 
(Fig. I12p . showing a power law behaviour with similar ex- 
ponents. Again, M4 shows a higher exponent (0.82) due to 
its additional thrust support in the nose cone. There might 
also be a slightly faster propagation for the M3 jet, where 
the magnetic field is not too much below equipartition at 
the jet inlet (see Fig. [2]), although this might also be just a 
temporal effect due to the jet-cocoon vortex interaction. 

The cocoon width, in contrast, shows different power 
law exponents depending on the density contrast, after the 
start-up phase is over. We find exponents of 0.59 (Ml), 0.52 
(M2), 0.39 (M3) and 0.24 (M4L) for the different models 
(Table |3]). Thus there seems to be a clear trend of decreas- 
ing exponents for lower jet densities, which holds true for all 
our cocoon width measures (Fig. I13p . Widths approaching 
an asymptotic value might mimic a similar behaviour, but 
so far, this is beyond our simulation data (except for M4). 
It seems reasonable that this is due to less overpressured 
cocoons for lighter jets, as it is the c ocoon pressure that 
drive s the lateral cocoon expan sion l|Kaiser fc Alexandeij 
1 19971 : ICarvalho fc O'Deal l2002bl ). If the cocoon pressure 
equals the ambient pressure, the sideways expansion of the 
cocoon would come to an end. 

Another consequence of this is the lateral bow shocks, 
compared to the corresponding cocoon width, being much 
further away for light jets (FiE. ll4|) . as found bv lZanni et al.l 
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Figure 12. Time evolution of the cocoon: full length and average 
width as function of time. MHD models in thick lines, correspond- 
ing hydro models in thin lines. 
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Figure 13. Power law exponents for cocoon widths as function 
of density contrast. Fits considered only data points where width 
^ 5 kpc. 



l|2003l ). too. Hence, except for the Hl/Ml models, the thick 
layer of shocked ambient gas grows continuously. 

The expansion of the cocoon for M4 is much different. 
After the initial phase the cocoon width settles down to a 
constant value and does not grow anymore. This is a conse- 
quence of the suppressed backflow in the nose cone, which 
then cannot inflate the cocoon anymore. 

All simulations with non-dominant magnetic fields show 
pronounced turbulence in their cocoons. This is evident from 
Fig. 1151 which shows the vector fields of velocity and poloidal 
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Figure 14. Bow shock / cocoon width ratio over time for the dif- 
ferent simulations as function of the full jet length. MHD models 
in thick lines, corresponding hydro models in thin lines. 



magnetic fields in LIC (line integral convolution) re presen- 
tation. The LIC technique jCabral fc Leedo"mlll993l ) allows 
the fine-grained depiction of vector fields, especially suit- 
able for turbulence, where structures even on smallest scales 
are present due to the turbulent cascade. We extended this 
to show the field magnitude, additionally, decomposing the 
information into brightness (showing the field direction as 
stream lines) and colour (field magnitude) in IfLS colour 
space. 

Cocoo n turbulence is drive n by quasi-periodic "vortex 
shedding" (jNorman et al.lll982l ) in the jet head, which in- 
jects vortices into the the cocoon. As these vortices move 
around and interact, vortex shedding affects the whole co- 
coon and drives its turbulence. While it occurs in our heav- 
ier jets, too, narrow cocoons suppress vortex interaction and 
the establishment of turbulence. We note that there may be 
feedback on the driving mechanism, as cocoon vortices per- 
turb the jet beam and thus influence the vortex shedding 
process itself. 

3.3.1 Aspect Ratio 

A characteristic property of the bow shock or cocoon is their 
aspect ratio TZ = length/width (Fig. [16}. Dependent on the 
density contrast, after a short initial phase of spherical ex- 
pansion (TZ ~ 1), the bow shock aspect ratios grow but 
converge for large bow shock diameters, approaching 1 for 
lighter jets (7^ = 1.4 for M3 and 1.1 for M4L). This means, 
the bow shock approaches a spherical shape for very fight 
jets. Once again, M4 is different, as the propagation in axial 
direction is faster, yielding significantly higher aspect ratios 
then M4L. 

The aspect ratios for the cocoons generally increase 
with jet length and are at early times systematically lower 
for the lighter jets. However, the light jets soon increase their 
aspect ratio (earlier for lighter jets) and then at later times, 
show aspect ratios even higher than their heavy counter- 
parts. As for the cocoon width evolution, we argue that this 
is due to cocoons, which come to pressure balance with the 
ambient gas earlier, so that lateral cocoon expansion stalls, 
but the axial propagation is still growing self-similarly. By 
comparing Fig. [8] and Fig. [16] one sees that once a source 
approaches pressure balance with its environment, it drops 
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Figure 15. Velocity field (upper panel) and poloidal magnetic 
field (lower panel) around the jet head, displayed in LIC (line 
integral convolution) representation to show the small-scale vec- 
tor field structure. The colours show the vector field magnitude, 
while the brightness modulation shows the field lines. 
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Figure 16. Aspect ratios TZ of bow shock (upper panel) and 
cocoon (lower panel) over full length. For the cocoon aspect ratios, 
the hydro models are omitted for clarity. 



out of self-similarity and increases its cocoon aspect ratio. 
For M3 this happens already early, while it does not happen 
for Ml until the end of the simulation. 

3.4 Entrainment 

The jet backflow at the contact surface between the cocoon 
and the ambient gas makes it Kelvin-Helmholtz unstable, 
and thus creates fingers of dense ambient matter that reach 
into the cocoon and are entrained. In numerical simulations 
this entrained gas is additionally mixed with the jet plasma 
due to finite numerical resolution. The amount of entrain- 
ment can be measured in terms of the cocoon mass since the 
mass of jet plasma usually is small compared to the mea- 
sured cocoon mass. Although the exact numbers depend on 
the exact cocoon measurement definition (Sect. I3.2|l . this 
seems to be a reasonably robust method. 

Figure [17] shows the time evolution of the cocoon mass. 
The entrained mass grows with a power law exponent only 
slightly below the exponent of the cocoon volume, showing 
a slowly decreasing but roughly constant fraction (5-10 per 
cent) of the initial mass in the occupied volume. However, 
there is no difference visible between purely hydrodynamic 
and MHD simulations in the entrained mass, as would have 
been expected. The reason for this is the missing stabiliza- 
tion of the contact surface, which is discussed later. However, 
it is evident from M3 in Fig. [T8]that the entrainment in the 
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Figure 17. Evolution of the cocoon mass as measure for the 
entrainment of dense ambient gas. The long-dashed line gives the 
injected jet mass Mj for the M2 run for comparison (Mj oc rjt). 
MHD models as thick lines, corresponding hydro models as thin 
lines. 



jet head is significantly smaller: the mass in a cylindrical vol- 
ume {Z € [—45, —35] kpc, radius 4 kpc) in the head region 
of 1VI3 is 3 X 10^ Mq without magnetic fields (H3) , compared 
to 9 X 10* Mq in the magnetized case, which is more than 
a factor of 3 lower. Hence entrainment is significantly sup- 
pressed in the jet head, but no change could be measured 
regarding the whole cocoon volume. 
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Figure 18. Comparison of liydro (H3) and MHD (M3) simula- 
tions at 15 Myr. The jet liead region is much more pronounced, 
Kelvin— Helmholtz instabilities are damped, and the jet has prop- 
agated a bit further. 
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Figure 19. Gains in thermal energy as fraction of the measured 
total injected energy, as function of the full jet length. 



3.5 Energy budget 

From the quick balancing of pressure within the cocoon, one 
might expect a strong conversion of (kinetic) jet power to 
thermal energy. This, in fact, is measured for our simula- 
tions. 

Figure [T5] shows the increase in thermal energy as frac- 
tion of the total injected power. Already for the heaviest jet 
(Ml), most of the injected (kinetic) power appears as ther- 
mal energy due to compression and irreversible entropy gen- 
eration at shocks. The thermal fraction increases not only 
with time, but is also much stronger for the lighter jets, 
where a thermalization of > 80 per cent is reached. Half of 
the thermal energy gain is found in the c o coon and half in 
the (shocked) ambient gas. lO'NeiU et all l|2005h find ^ 40 
per cent of the jet power in the thermal ambient gas for 
their 3D jets with density contrast 0.01 in a uniform atmo- 
sphere, while in our simulations we find ~ 35 per cent, which 
is quite good agreement. 

Magnetic energy only has a very small contribution (be- 
low 1 per cent), except for M4, which is magnetically dom- 
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Figure 20. Evolution of the magnetic energy with time. The 
long-dashed line shows the injected amount of magnetic energy. 
For the case of M4L, magnetic energy was multiplied by 100 to 
account for the 10 times smaller field strengths. 



inated and which has a magnetic energy contribution ris- 
ing up to 5 per cent. More than 90 per cent of the mag- 
netic energy is located in the cocoon. For all runs except 
M4, the magnetic energy that is actually measured is signif- 
icantly larger than the injected magnetic energy (this effect 
is stronger for the lighter jets) and it grows faster then just 
linearly in time - approximately with a power law expo- 
nent of 1.2 (Fig. I20|) . Hence, other forms of energy seem to 
be converted into magnetic energy. For M4, the measured 
magnetic energy is lower than expected from the nozzle val- 
ues, which may indicate that the additional thrust in the 
nose cone actually consumes magnetic energy. 

The remaining fraction is kinetic energy, which is de- 
creasing more and more with lower jet density. Although 
the jet beam is the only energy input to the system (at the 
nozzle) , its contribution to the total kinetic energy is 10 per 
cent or less, and 50 to 30 per cent is in the cocoon. The re- 
mainder, 50 to 70 per cent, comes from the outward moving 
shocked ambient gas. 



3.6 Magnetic fields 

Magnetic fields are not only passive properties of the jet 
plasma, but an active ingredient for the dynamics. One pa- 
rameter describing this is the ratio between thermal and 
magnetic pressure of the plasma {(3 = Snp/B'^). For the sim- 
ulations described here, we used a fixed value for jet speed, 
Mach number and magnetic field. Thus, the plasma /3 can- 
not be constant throughout the different runs (see Tab. [2]). 
While Ml and M2 have passive magnetic fields, M3 and M4L 
have fields with significant impact, and for M4 they are even 
dominant. 

The helical field configuration in the jet initiates an in- 
triguing interplay between kinetic and magnetic energy. Al- 
though the jet matter is injected without any rotation, the 
Lorentz force from the helical field generates a toroid al ve- 
locity component, as also found bv lKossl et al.l ()l990l ). This 
effect is stronger for the runs with stronger magnetic fields 
(lower plasma P). The rotation does not originate from per- 
sisting angular momentum from the jet formation, which 
should be very small due to the expansion of the jet. Also, 
it is not continuous throughout the beam and even chang- 
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ing sign at some internal shocks and interaction with cocoon 
vortices. 

When the plasma reaches the terminal shock, it flows 
away from the axis radially and turns back, forming the 
backflow that inflates the cocoon. Rough conservation of 
angular momentum I then produces a radially declining an- 
gular velocity fl = l/B? (differential rotation). Writing the 
induction equation in cylindrical coordinates, 

^ = -i? (^p ■ V) - B^V ■Up + R{Bp- V) n, (2) 

it becomes evident that this shearing transforms poloidal 
field Bp into toroidal field B^, also transferring kinetic en- 
ergy into magnetic energy, which explains why the contri- 
bution of magnetic fields to the total energy is higher than 
its injected contribution. 

We note that this creation of toroidal field in the jet 
head is not an artefact of axisymmetry, but merely a conse- 
quence of allowing threedimensional vectors in the simula- 
tion {u and B). We do not expect this to be much different 
in full 3D, apart from a naturally more complex structure in 
the details. What, in contrast, most probably is an artefact 
of axisymmetry is the persistence of the toroidal field com- 
ponent in the cocoon. The cocoon plasma is highly turbulent 
(Sect. I3.3.3|) with relatively little systematic motion, which 
is an intrinsically threedimensional phenomenon. This can 
easily convert toroidal and poloidal field into one another, es- 
tabhshing some dynamical equilibrium between those com- 
ponents, but maintaining the overall field strength. 

Comparing purely hydrodynamical models with the 
MHD models (Fig. llSp . we find that the global properties, 
such as bow shock and cocoon sizes, are generally robust if 
the magnetic fields are not dominant (as with M4) . The de- 
tails, though, are different. While the hydro models show a 
ragged contact surface between jet plasma and ambient gas 
due to Kelvin-Helmholtz (KH) instabilities excited by the 
backfiow, the MHD runs show a pronounced jet head, which 
is clearly more stable, si nce the KH instability is damped by 
the mag netic fields (e.g. lMiura fc Pritchett ''l982^. Magnetic 
tension acts as a restoring force on the growing instabilities, 
suppressing entrainment of ambient matter and "fingers" of 
dense gas reaching into the backfiow, which is evident from 
the clearly lower average density in the jet head region. The 
stabilizing effect appears at /3 ~ 10 in the jet head. For the 
simulations with weaker fields there is no noticeable differ- 
ence between the magnetized and the pure hydrodynamics 
case. 

Damping of the KH instability by magnetic fields, how- 
ever, only works with the field component parallel to the 
instability wave vector, which in turn means that in ax- 
isymmetry only the poloidal magnetic field can damp the 
instabilities at the contact surface. 

Although the earlier-mentioned shearing mechanism 
amplifies magnetic fields and should therefore provide even 
more damping of KH instabilities, we cannot see this effect 
further away from the jet head, because in axisymmetry the 
backward reaction (toroidal to poloidal) cannot work and 
thus the poloidal component becomes too weak (Fig. I15|) . 
As the magnitude of the magnetic field in the cocoon is 
as strong as in the jet head, it seems reasonable that with 
balanced magnetic field components in reality, the contact 
surface could be stabilized. 



The toroidal field B^ is directly related to the generat- 
ing current jp, which is shown in Fig. [26] as field lines. Our 
toroidal field setup describes a situation where the poloidal 
currents leave the nozzle axially in the jet core, turning back 
in the sheath. As the backflow develops, the poloidal current 
flows along the contact surface with typical integrated cur- 
rents of several 10^* amperes (|Camenzindlll990l : BlandfordI 
[2008!). The toroidal field in the cocoon, built-up by the shear- 
ing in the jet head, seems to form its own current circuits. 
The gross radial behaviour cx R (Fig. I22|l can be at- 
tributed to the relatively uniform distribution of the axial 
current through the planes perpendicular to the jet beam. 

If the toroidal field is strong in the jet head region, the 
Lorentz force /l — jxB/c produces additional thrust for the 
jet propagation due to the strong radial current component, 
which is evident for M4, showing a pronounced nose cone, 
and may also explain the slightly faster propagation of M3 
with respect to H3 (Fig. [T2)) . 

Inside the beam, the magnetic field stays mostly 
poloidal, as injected, but near the terminal shock it is com- 
pressed axially, directed off the axis and sheared, producing 
strong toroidal field loops (Fig. I27|l . 

Finally, we turn to volume-weighted 2D histograms of 
magnetic pressure pm = B^/Stv and thermal gas pressure p 
in Fig. 1281 where the contributions from only the jet beam 
and all the jet plasma is shown separately. 

The jet nozzle is located at (logp, logpm) ~ 
( — 10, —10.5) as a vertical line (constant pressure, but ra- 
dially varying magnetic field). As the matter fiows through 
the beam, internal shocks (cf. Figs. [6land ll8|) cause strong 
changes in pressure whereas the plasma /3 remains un- 
changed (magnetic field is compressed with the plasma), 
leading to lines originating from the nozzle location parallel 
to the overplotted /3 — const lines. The plasma /3 somewhat 
increases along the beam when it it interacts with the co- 
coon vortices, and thus creates some down-shifted parallels. 
There is no clear separation between the beam and the en- 
closing cocoon in the beam-only diagram, hence both shear 
layers of the beam and cocoon gas are contained in the wide 
area below P = 10. Still, there are strong pressure changes 
indicated by the wide horizontal distribution. 

The distribution of cocoon cells widely spreads both 
to higher and lower magnetic fields from this area. The pro- 
nounced trail downwards is the transition to the ambient gas 
through entrainment; since the ambient gas is essentially un- 
magnetized, it is located even below the lower border of the 
figure (Fig. I28p. The radial increase of magnetic field in the 
cocoon yields the extension towards lower plasma /3 (see also 
Fig. [25}. The spiky features around /3 ~ 2 are single vortices 
in the outer parts of the cocoon, where the pressure drops 
towards the centre due to centrifugal forces together with 
a slight increase of toroidal field. Altogether, the spread of 
the cocoon cells is considerably larger in magnetic pressure 
than in thermal pressure. 

The situation shown in Figs. [28] and [29] is typical for 
the time evolution of these diagrams. Clearly, some fea- 
tures are appearing, changing and disappearing continu- 
ously, such as individual internal shock lines or the cocoon 
vortex spikes. The general structures in the diagrams per- 
sist at all times. There are, however, two systematic changes 
with time: Firstly, the "cocoon bump" in Fig. [29] (log /9 ~ 1) 
grows due to cocoon expansion, eroding the "ambient bump" 
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Figure 21. Toroidal magnetic field magnitude for M3 at 15 Myr. The right jet has positive toroidal field, the left jet has negative sign. 
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Figure 22. Toroidal field in three slices parallel to the mid- 
plane, and the average of the whole region Z = 10 ... 30 kpc. 
Outside the jet radius rj = 1 kpc, the toroidal field shows a 
roughly linear increase with R. Although in every slice drops 
to at -R 1 kpc, the average does not drop to (and thus is 
a bit misleading), as the positions are slightly offset for different 
slices. 

(log/3 ~ 9), and moves to the left, faster at early times and 
then becoming continuously slower. Secondly, as the cocoon 
pressure drops, the cocoon distribution of Fig. l28l moves to- 
wards the left (and somewhat down due to the mostly con- 
stant distribution in /3 at late times) , and grows with cocoon 
volume, too. 



3.7 The lightest jets 

The lightest jet in the series, M4, shows a very different be- 
haviour from the other runs due to its strong magnetic fields, 
thus a run with lower magnetic fields (M4L) was performed 
in addition. In this subsection, we focus on the specific prop- 
erties of and differences between these two runs. 

Both simulations show unstable beams, which are tem- 
porally stopped, deflected or disrupted. This is quite natural 
for the very light jets, where the impact of cocoon vortices 
hitting the beam is stronger, when the beam shows lower 
inertia but the cocoon gas is dense due to entrainment and 
mixing with the dense ambient gas. This destabilization is 
particularly strong in axisymmetry, as the vortices cannot 
"miss" the beam as they could in 3D. For M4L, after a 



strong deflection of the right beam near the nozzle {t ~ 20 
Myr), a small region with strong poloidal held piles up just 
next to the nozzle and creates a magnetic layer {/3 > 1) 
at the beam boundary. This protects the beam from cocoon 
vortices and entrainment, and from there on inhibits disrup- 
tions of the right jet, which then propagates more quickly 
than the left jet. At the end of the simulation, the right jet 
is ~ 20 per cent longer than the left jet and shows an al- 
most undisturbed beam up to the jet head. More detailed 
examination of this phenomenon may be interesting, but as 
it was orfly introduced by chance, the details are difficult 
to reproduce and beyond the scope of this paper. None the 
less, the overall propagation of the jet within the simulated 
time (Sect. [3T3.3|I is not much affected by this. 

Keeping the jet speed and the Mach number fixed, the 
ratio of the thermal pressures of ambient gas and jet noz- 
zle changes with density contrast, yielding an underpres- 
sured jet for M4 and M4L. For M4, the magnetic field in 
the nozzle is already stronger than equipartition and the 
Alfven speed is higher than the sound speed. This run is 
dynamically dominated by the magnetic field and shows a 
pronounced nose c one, which is known for jets with strong 
toroidal fields (Cl arke et al.lll98^ ). Magnetic tension pinches 
the jet matter into an narrow tube of 2.5 to 3.5 kpc ra- 
dius, completely suppressing a backfiow and thus prevent- 
ing the formation of a wide cocoon. The simple case of a 
plasma column in radial magnetostatic equilibrium keeps 
p + B^/A-K constant. If B'^/A-k approaches the thermal pres- 
sure, the magnetic pinch becomes important. In our case, 
the toroidal field in the plasma column is relatively homo- 
geneous, showing a (volume-weighted) distribution mostly 
between 30 and 50 \iG, while the thermal pressure lies (ra- 
dially decreasing) in the range 1 ... 3 x 10"^" dyne/cm^, thus 
matching p ~ B'^/Att and being just around equipartition. 
These values are not the ones set by the jet nozzle, although 
those obey p ~ B'^/Atv, too. The twisting and shearing pro- 
cesses described in the previous subsection are very strong 
due to the equipartition-level magnetic fields, the rotation 
around the jet axis can make up a large fraction of the to- 
tal velocity, and the toroidal field component grows to the 
measured values. 

iKrause fc Camenzindl l|200ll) examined the convergence 
of a nose cone simulation and found that the Mach disc 
retreated towards the nozzle and thus did not converge. Also 
in M4, the Mach disc is very near to the jet nozzle, and 
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Figure 23. Fractional contribution of toroidal field to the total magnetic field. Simulation M3 at 15 Myr. 
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Figure 24. Angular velocity |n| = \v^\/R in units of s~^, scaled 
logarithmically. Note that there is a strong decline from the beam 
to the backfiow. 



the velocities after that shock are subsonic (although the 
nose cone itself propagates faster than the jet head in M4L). 
Thus, it is unclear, how reliable the run is. We also n ote that 
the m agnetic pinch is subject to MHD instabilities (jClarkd 
1 19931 ). which might produce blobs and disrupt the plasma 
column in 3D. However, as this nose cone is produced by 
the magnetic tension of the strong toroidal field, this is not 
applicable to strong poloidal fields, which cannot provide the 
necessary hoop stress, although it seems difficult to maintain 
a strong poloidal field along an interacting beam without 
converting part of it into toroidal field, which then might 
again pinch the plasma. 



4 DISCUSSION 
4.1 Magnetic Fields 

Effects of magnetic fields naturally depend on their strength. 
Trying to understand the smoothness of jet cocoons in 
galaxy clusters, we concentrated on magnetic fields in jets 
which are not dominant, but still have significant effects 
on the jet dynamics, the best example for this being the 
M3 run with average p lasma l3 — 8. It is well known 
ijMiura fc Pritchettlll982l ) that magnetic tension can damp 
or suppress Kelvin-Helmholtz (KH) instabilities and hence 
it may be the key to stabilizing the contact discontinuity be- 



tween jet and ambient gas. However, how this applies to the 
complex case of jet-ambient interaction is not yet known. 

We emphasize that much care was taken to use a glob- 
ally consistent setup for our very light jets, in particular: 
keeping the bow shock inside the computational domain at 
all times; simulating bipolar (back-to-back) jets to remove 
an artificial boundary condition in the midplane and allow 
interaction of the backflows for a realistic lateral expansion; 
and using a configuration which confines the magnetic field 
to the jet and has closed field lines instead of a homogeneous 
magnetic field reaching to infinity, which is then effectively 
anchored in the ambient gas. The assumed simplifications, 
axisymmetry and a constant ambient density, make extrac- 
tion of the underlying physics easier and effects of relax- 
ing tlwseforh^^rodynamic jets were previously investigated 
bv lKraus3 l|2005l ). Thus we expect to at least qualitatively 
model the situation realistically. 

Two main effects arise from the inclusion of magnetic 
fields: Firstly, in the jet head, we see that the provided mag- 
netic fields in the jet do indeed stabilize the contact surface, 
which produces a pronounced jet head and lobe s, similar to 
the ones seen in Cygnus A (|Carilli et al.l Il99lh and other 
classical double radio sources. Effects from an ambient den- 
sity profile can be excluded due to the prescribed constant 
density atmosphere. Furthermore, the entrainment of ambi- 
ent gas is significantly smaller there than without magnetic 
fields. 

Secondly, jets prove to be efficient generators of mag- 
netic energy, transferring part of their huge kinetic power to 
magnetic fields through shearing in the jet head. This relies 
on some rotation of the beam plasma, which will (as seen in 
the simulations) generally be present for a non-zero toroidal 
field component. Some t oroidal field is expected if the mostly 
axial field in the beam l|Bridle fc Perlevlll984) is perturbed 
three-dimensionally and from jet formation models, where 
the toroidal field is necessary for jet coUimation at least at 
small scales. The shearing mechanism provides a source of 
magnetic energy for the cocoon and furthermore affects the 
magnetic field structure at the hotspots and possibly some 
internal shocks. A radial and toroidal field component in the 
beam is known to be compressed by the terminal shock and 
is then visible as a strong magnetic field perpendicular to 
the jet axis. The jet head shearing provides another mecha- 
nism, independent from compression, to greatly enhance the 
toroidal field and thus produce a perpendicular field com- 
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Figure 26. Poloidal current field fines (contours of RB^). Grey: negative values, black: positive values. Model M3 at t = 15 Myr. 



ponent stronger than that expected from compression. For 
jets pointing more towards the observer, the toroidal field 
around the hotspot region may become observable. 

This may be relevant for several observational findings, 
one being the smoothness of radio cocoons. We have shown 
that even if the plasma /3 is of order ten, only, the fields 
in the backfiow and the cocoon respectively will be strong 
enough to damp KH instabilities at the cocoon-ambient gas 
interface and yield a morphology much smoother than seen 
in hydrodynamic simulations, recon ciling simulations with 
observations of sources as Cygnus A jLazio et al.ll2006l). Pic- 
tor A jPerlev et al.lll997h or Hercules A l|Gizani fc Leahvl 
I2OO3I ). where the latter seems to be a past high-power 
source. Due to the 2.5D nature of the simulations, the ef- 
fect is restricted to the jet head region. In a full 3D sim- 
ulation, we expect therefore the cocoon-ambient interface 
to be more stable even further back from the hotspots. 
The amplification of beam magnetic fields in the "jet head 
machine" furthermore is consistent with the observation of 
mag netic fields in the cocoon just s omewhat below equi parti- 
tion (jHardcastle fc CrostoD|[2005l : iMigliori et al.ll2007l ). Ad- 
ditionally, the magnetic field predominantly perpendicular 
to the jet axis in weak FR I sources might be related to 
the expansion of the jet, which by the shearing would create 
strong toroidal fields in the absence of strong turbulence. 
Even though the beam rotation can change much due to in- 
teraction with the cocoon and shocks and even change sign. 



the helicity of the toroidal field is not changed and can thus 
link the fiel d at large scales with the field topology near the 
black hole ijCabuzda et al.ll2008l ). 

For the magnetic field topology in the rest of the co- 
coon, axisymmetry is a major limitation, contrary to the 
effects discussed before. Magnetic field in a toroidal con- 
figuration cannot damp KH instabilities in axisymmetry 
since no magnetic tension is available as restoring force, 
while poloidal field could do so. Fortunately, the jet head- 
generated toroidal field in the turbulent cocoon partly would 
be converted into poloidal field in three dimensions, estab- 
lishing some dynamical equilibrium between the components 
but keeping the overall field magnitude or amplify it even 
further, and this makes the cocoon magnetic field a rea- 
sonable explanation for the smooth contact surfaces. As a 
future step we will examine this effect in three dimensions to 
be able to quantify the amount of damping and suppressed 
entrainment of ambient gas in the cocoon. 

However, despite the inability to actually produce the 
expected smooth contact surfaces in axisymmetry away from 
the jet head region, there is no reason to assume that the am- 
plification of magnetic fields should be in three dimensions 
any different than shown in our simulations, as the plasma 
dynamics is not very different and the shearing mechanism 
in the jet head simply relies on the off-axis flow of plasma, 
which also happens in 3D. Furthermore, we are not aware of 
any reason that the fleld magnitude in the cocoon should be 
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Figure 27. 3D magnetic field lines from within the beam with a transparent grey-scale logarithmic density slice in a plane through the 
jet axis. Model M2 at t = 1.3 Myr. 



much different in 3D. It is unclear, though, how the spatial 
distribution of the magnetic field would look like: 3D tur- 
bulence might want to distribute field strength rather uni- 
formly in the cocoon, but formation of a large-scale poloidal 
current may try to establish a radially increasing toroidal 
field. Observations indicate that magnetic field strengths 
withi n the cocoon may vary considerably l|Goodger et al.l 
l2008l ). 

It may be interesting to note that the amplification of 
magnetic field is quite related to dynamo action as in the 
sun. The shearing {Q effect) is just the same and solar con- 
vection is replaced by jet-driven cocoon turbulence, but the 
location of these actions are different and they are exter- 
nally powered (by the beam thrust) instead of self-sustained. 
The spatial separation of the two effects and the (at least 
roughly) isotropic turbulence, however, prevent the forma- 
tion of an outstanding large-scale poloidal field. 

The uncertainty in the magnetic field topology in the 
cocoon also applies to the distribution of plasma P in the 
system. We (expectedly) found that plasma (3 is unchanged 
throughout shocks despite a gradual increase along the beam 
(which might also be due to limited resolution of the beam 
and entrainment). Thus the assumption of a fixed fraction 
of equipartition to generate synchrotron emission maps from 
hydro simulations seems to be quite justified. However, this 
was not found to be true for the cocoon, where a wide dis- 
tribution of P was found and deriving synchrotron emission 
from hydro models thus may be far away from MHD results. 



But as mentioned, this result is expected to change in 3D, 
apart from having relatively low /3 in the cocoon. Emission 
maps of our simulations and comparison to hydro models 
are beyond the scope of this work and will be presented in 
a subsequent paper. 

The amplification of magnetic fields is also particularly 
intere sting for the que stion of the origin of lobe magnetic 
fields. iDe Youn3 (|2002l ) pointed out that equipartion fields 
in the lobes cannot be passively advected with the plasma 
from the jet beam due to fiux conservation arguments. The 
beam magnetic fields would have to be of order 0.01 G or 
higher, certainly above equipartion, which would result in 
enormous synchrotron losses, luminosities incompatible with 
observational limits and probable disruption of the jet due 
to the magnetic pressure. Hence, the magnetic field must 
be amplified by some mechanism, and lDe Yound argues for 
turbulent amplification in the hotspot fiow, though it is not 
easy to meet the necessary requirements for this. The shear- 
ing in the jet head, which is seen in our simulations, in con- 
trast, almost inevitably provides this amplification and can 
therefore explain the strong lobe magnetic fields or at least 
contribute to their field strength. In fact, the simulations ex- 
hibit field magnitudes in the cocoon that are comparable to 
field magnitudes in the beam and consequently have similar 
plasma /3 since the beam and cocoon pressures came to bal- 
ance. We conclude that shearing due to off-axis fiow of the 
plasma provides a natural explanation for the lobe magnetic 
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Figure 28. Pm-p histogram for tlie beam (defined by kinetic 
energy flux ^ f % of the maximum value) and all the jet matter 
from simulation M3 after 15 Myr. Lines of constant plasma /3 
with values of 1, 10, 100, 1000 and 10 000 are overlaid. 
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Figure 29. Volume-weighted histogram of plasma fi for M3 at 
0.5, 7.5 and 15 Myr. 



fields and allows equipartition jets to inflate an equipartition 
cocoon. 



4.2 Dynamical Evolution 

X-ray observations of the ambient cluster gas contain valu- 
able information about several jet and AGN properties and 
self-similar models can give easy access to underlying phys- 
ical parameters. In the present paper, we are able to con- 
firm agreement of our numerical simulations with self-similar 



models llFalld Il99ll: iBegelmaiil 1 19961 : iKaiser fc Alexandeij 
1 19971 : iKomissarov fc Falld Il998l ) for the bow shock propa- 
gation. Excentricity of the bow shock and its Mach number 
seem to be an easy way to compare theoretical models with 
observations, without the need for uncertain assumptions on 
the emission of the radio plasma. 

The weak and roundish bow shocks in observations in- 
dicate that models of very hght jets (with density ratios 
< 10~^) are necessary for most cluster sources. Although 
we chose a simplified setup with a constant ambient gas and 
axisymmetry, the simulations are in the regime of observed 
values for various sources and self-similar models generalize 
this behaviour for declining cluster profiles, which was al- 
ready examined for very hght hydrodynamic jets bv lKrausd 
l|2005h . As our runs, with the clear exception of the magneti- 
caUy dominated M4, propagate as their hydro counterparts, 
only minor deviations from those results are expected, ex- 
cept where specific source properties are to be included. 

Contrary to the bow shocks and the jet length, we 
find that the cocoon width in general does not evolve self- 
similarly but for lighter jets grows with lower power law 
exponents and the mean cocoon pressure drops more slowly 
than expected. A lthough this may seem une xpected, it was 
already stated bv lKaiser fc Alexandeij (119971 ), that contrary 
to the bow shock, the self-similar evolution of the cocoon 
depends on the physical model for the post-hotspot fiow 
and thus, deviations are to be expected if these assumption 
do not hold in the simulations. Since very hght jet cocoons 
are less overpressured and approach the ambient pressure 
sooner, the sideways expansion becomes slower and may 
even stall, letting their aspect rat io (length to width ratio) 
grow. This is in fact observed bv lMuUin et al.l l|2008l ). who 
find a wider range of aspect ratios, once the source size ap- 
proaches ~ 100 kpc. Similar behaviour would be expected 
for the heavy jets, although at much later times. Thus, co- 
coon evolution depends sensitively on the question of over- 
pressure, which can be addressed by the strength of the 
lateral bow shock. Self-similar models, in contrast, assume 
that the ambient pressure is neglig ible. H omissarov fc Falld 
(1998) defined two scales, 1^ and Lc, between which they 
found self-similar evolution. The lower bound Zc, where the 
swept-up mass equals the jet mass, is much smaller than 
the jet radius in our simulations, and the upper bound Lc, 
where the ambient pressure becomes important, is compa- 
rable to our computational domain size. Accordingly, while 
they observe the self-similarity being established, we observe 
its end, explaining why our less overpressured numerical so- 
lutions gradually deviate from a self-similar evolution. 

Furthermore, cluster density profiles make cylindrical 
cocoons rather than elliptical on es due to the weaker den- 
sity contrast at larger distances (lKrausdl2005l ). Altogether, 
this makes us confident that our simulations reasonably well 
describe observed cluster sources. 

In contrast to bow shocks, measurements of the co- 
coon shape are complicated by cooling of the relativistic 
electrons, which limits observations to the outermost parts 
(lobes). While radio observations show the high-energy elec- 
trons in the cocoon as lobes, single-fluid MHD simulations 
only trace the low-frequency emitting matter and can only 
show the low-freque ncy radio morphol ogies (cf. high and low 
frequency images in lCarilli et al.|[l99ll ), which generally suf- 
fer from low spatial resolution. This situation fortunately 
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will much improve in the future with new telescopes as LO- 
FAR or the SKA, which will allow more detailed studies of 
cocoon dynamics and turbulence. Until then, X-ray images 
of cavities and (in some cases) the inverse-Compton emis- 
sion off the cosmic microwave background may supplement 
available low- frequency radio maps. 

Scheuer (1982) introduced the "dentist's drill" to refer 
to a moving working surface, which therefore widens the jet 
head and the lobes. Very light jets naturally show extensive 
cocoons and varying deflection of the beam widens the jet 
head and hence, even in axisymmetry, show something very 
similar to a "den tist's drill". While this does not exclude 
beam precession ijSteenbrugge fc Blundelll2008l ). it is does 
not require it and no large precession amplitudes are needed 
a priori. 

We expect for multiple outbursts of different power 
in the same cluster, indicated by "ghost cavities" (e.g. 
iFabian et alll2006l : IWise et al.1 l2007l 'l. that their evolution 
crucially depends on the history of the past outbursts, as 
these push the dense cluster gas aside, letting the new out- 
burst propagate with different density contrast. In this case, 
the new jet might quickly push forward to the old jet size, 
then resuming its work on the dense ambient gas. The mor- 
phology of the cavities may allow the determination of the 
respective density contrasts and thus could shed light on the 
outburst history. 

The thermal interaction of jets with the intra-cluster 
medium is less accessible to direct comparison with ob- 
servations. Slower jet head propagation is responsible for 
the strong impact of the beam at the working surface and 
a high thermalization; some conversion of kinetic to ther- 
mal energy will additionally occur near or in the beam 
due to beam destabilization, but may be less effective in 
3D. Despite the dominant power source is the kinetic jet 
power, this strong thermalization converts most of the in- 
put power to thermal energy - about half of this in the 
shocked ambient gas and half in a cocoon filled with high- 
entropy plasma, which eventually may transfer at least part 
of its energy to the entrained cluste r gas. This is in line 
with findings of o t her authors (e.g. [Reynolds et al.l |2002| : 
lO'Neill et ;il]|2005l : IZanni et al] 120051 ), where the latter au- 
thors conclude that up to 75 per cent of the energy can be 
dissipated irreversibly and thus is available for heating in the 
intra-cluster medium, as required by t he X-ray luminosity- 
temperature relation (|Magliocchetti fc BriiggenI 120071 ') and 
to pro vide "radio-mode" f eedback for models of galaxy evo- 
lution (|Croton et al.ll2006l ). 

Since only the hot gas phase is simulated, effects on 
the cold or warm phases of the interstellar medium (ISM) 
of galaxies are difficult to estimate. Clearly, the thermal- 
ization efficiency cannot be simply applied to the cold gas. 
Simulations of multi-phase turbulence in the jet cocoon by 
iKrause fc Alexand er (2007) with their higher spatial resolu- 
tion can resolve the different phases and provide a comple- 
mentary view ( "microphysics" ) onto the jet-cloud interac- 
tion. However, even if the thermal energy mostly is deposited 
on the hot gas phase (at larger distances), it is evident from 
our simulations that the jet cocoon is a rich reservoir of tur- 
bulent kinetic energy which will act on the cold gas phase of 
the galaxy over a time scale corresponding to the decay time 
scale of the cocoon turbulence. For a jet of power 10*® erg 
s~^ being active for 10^ years, the turbulent energy stored 



in the cocoon is expected to be on the order of a few x 10 
ergs and over a time possibly longer than the jet activity 
will interact with the cold ISM phases. 

Another interesting result of the present simulations is 
the excitation of sound waves in the ambient gas by vor- 
tices in the turbulent cocoon, which is more effective for the 
very light jets with their extended cocoons. Vortex shed- 
ding ([Norman et al . 1982) quasi-periodically occurs in the 
jet head, and the vortices then are advected with the back- 
ffow into the cocoon and provide an intermittent source for 
the turbulent cascade, producing pressur e waves. Waves like 
these are visible in the Per seus cluster (|Fabian et al. I I2OO6I : 
IShabala fc Alexandeij[2007l ) and, although being hard to ob- 
serve, may be an ubiquitous feature in galaxy clusters with 
current or past jet activity. Their typical wave length might 
yield a link to jet dynamics and cocoon turbulence. In the 
lightest of our jets (M4L), the bow shock is just about to 
turn into a sound wave and then simply would join the en- 
closed sound waves. Viscous damping may be a mechanism 
to reduce the amplitudes in addition to the growing wave 
area and is another candidate for preventing cooling ffows 
(Fabian et al. 2005), but in our scenario would be related to 
the jets rather than to the AGN itself. 

Axisymmetry naturally imposes some constraints on 
the dynamics, which have to be considered carefully. Jet 
beams in high-power sources are essentially axisymmetric 
objects and effects of the full third dimension are merely 
perturbations from axisymmetry. However, this obviously 
is not true when beam stability or non-axisymmetric effects 
are explored specifically. While generally 3D jets are subject 
to a greater number of instabilities, for very light jets there 
is an opposing effect of an increased number of dimensions. 
While in 3D, cocoon vortices often will miss the beam or are 
slightly deflected, this is not possible in axisymmetry and 
the beam thus is destabilized, deflected or disrupted more 
easily which is most evident from ou r light est run (M4L). As 
seen in the very light jets of lKraus 3(|2005D. the beam stabil- 
ity improves when going to full three dimensions. For most 
results, however, energetics and scaling behaviour are not 
expected to change signiflcantly in 3D, notable exceptions 
to this being cocoon turbulence, magnetic field topology and 
stability of the contact discontinuity. 

Cocoon turbulence further away from the jet head cer- 
tainly will differ with increased dimensionality as the in- 
creased number of degrees of freedom for vortices allow them 
to turn in all directions and interactions between colliding 
vortices will be different. Though, we expect the effects on 
cocoon morphology to be within reasonable limits, as the ki- 
netic energy in the cocoon is lower than the thermal energy 
by factors of > 3 for rj ^ lO"'^ and hence, effects of thermal 
pressure will dominate. 



5 CONCLUSIONS AND SUMMARY 

We performed a series of axisymmetric hydrodynamic and 
magnetohydrodynamic simulations of bipolar very under- 
dense jets in a constant density atmosphere. The magnetic 
field is mostly confined to the jet with a helical topology. 

(1) We find that the magnetic fields damp Kelvin- 
Helmholtz instabilities in the jet head and stabilize it. They 
produce smoother and more pronounced outer lobes already 
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with a plasma /3 ~ 10. The entrainment of ambient gas into 
the cocoon is considerably suppressed there. This morphol- 
ogy is more consistent with observations of powerful dou- 
ble radio sources than arc hydrodynamic simulations, which 
show a ragged cocoon boundary. 

(2) Magnetic fields are efficiently amplified in the jet 
head by shearing as the plasma streams off the jet axis. This 
originates from a rotation of the beam which we find to be 
a general result of a toroidal field component being present 
(yet not necessarily dominant) in the jet. The shearing con- 
verts part of the huge kinetic energy into magnetic energy 
and provides the cocoon with a magnetic field much stronger 
than expected from flux conservation, in some regions 
even approacliiiig equipartition. These findings are consisis- 
tent with recent observations of near-equipartition magnetic 
fields in cocoons derived from radio/inverse-Compton emis- 
sion observations. Already in our axisymmetric simulations 
the fields are in principle strong enough to stabilize the con- 
tact surface between the cocoon and the ambient gas all- 
over the cocoon and not only in the jet head. The necessary 
change in field topology would be a consequence expected 
of fully threedimensional turbulence in the cocoon. 

(3) The amplified magnetic field is mostly toroidal, re- 
sulting in a stronger contribution of the field component 
perpendicular to the jet axis than expected from pure com- 
pression of magnetic fields at the hotspots. It is also ex- 
pected at locations where a jet widens considerably (as in 
FR I sources) . In the backflow and the cocoon, however, tur- 
bulence will probably establish some balance between the 
magnetic field components, which could not be established 
in axisymmetry. 

(4) The very light jets show round bow shocks with 
low Mach rmrnbcrs. We find that the bow shocks evolve 
self-similarly and hence give a simple link between obser- 
vations and some underlying physical parameters. The co- 
coon width, however, evolves self-similarly only for jets in 
their highly overpressured phases, but grows slower as the 
cocoons approach pressure balance with the ambient gas and 
the bow shock Mach number drops. These sources thus are 
surrounded by thick layers of shocked ambient gas. 

(5) The jet cocoon shows highly turbulent motion. It is 
driven by vortices shed in the jet head, which are advected 
with the backflow. Interaction of these vortices with the am- 
bient gas excites waves and ripples in the shocked ambient 
gas, which are joined by the dissolving bow shock at later 
stages. 

(6) The strong thermalization that occurs for very light 
jets transfers most of the jet power to the thermal energy of 
the cocoon and the shocked ambient gas, maJiing it available 
for heating of the cluster gas and radio-mode feedback. In 
addition to this, the turbulent motion in the cocoon is asso- 
ciated with a considerable amount of kinetic energy (~ 10 
per cent of the jet power) that may provide efficient feedback 
onto the cold phase of the galaxy's interstellar medium. 
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